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We report exact model calculations of the spin-dependent tunneling in double magnetic tunnel 
junctions in the presence of impurities in the well. We show that the impurity can tune selectively 
the spin channels giving rise to a wide variety of interesting and novel transport phenomena. The 
tunneling magnetoresistance, the spin polarization and the local current can be dramatically en- 
hanced or suppressed by impurities. The underlying mechanism is the impurity-induced shift of the 
quantum well states (QWS) which depends on the impurity potential, impurity position and the 
symmetry of the QWS. 
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Tunneling of spin-polarized electrons through mag- 
netic tunnel junctions (MTJ) has attracted0 wide and 
sustained interest in the past few years, both experimen- 
tally Ai^iS] and theoretically[30,|l,iS|- This is due 
to the potential applications of MTJ in spin-electronic 
devices, such as magnetic sensors and magnetic random- 
access memories. The key point for these applications 
is the tunnel magnetoresistance (TMR), i.e. the depen- 
dence of the tunneling current on the relative orientation 
of the magnetization of the ferromagnetic layers, which 
can be changed by an applied magnetic fieldQ- 

Double magnetic tunnel junctions (DMTJ) consist of 
a central metallic layer (quantum well) between two in- 
sulating barriers and two ferromagnetic electrodes. The 
insulating layers are thin enough for electrons to tun- 
nel through the barriers if a bias voltage is applied be- 
tween the electrodes. The TMR behavior in DMTJ is 
determined by quantum well states (QWS) formed in the 
middle layer when a resonance condition is fulfilled 0. 
The TMR can be dramatically enhanced when spin- 
polarized electrons resonantly tunnel through the mid- 
dle layer 0. Theoretical formulations of the TMR in 
DMTJ are usually based on models which assume per- 
fect systems 7, 9, 10, 11, 12, 13, 14]. The TMR exhibits 
an amplitude-varying oscillatory behavior as a function 
of the thickness of the middle layer with a period of 
ir/kp, where kp is the spin-dependent Fermi wave vec- 
tor in the middle layer ^^l- Chshiev et al introduced 'o'l 
phenomenologically the electron scattering in the mid- 
dle layer without taking into account vertex corrections. 
Consequently, they introduced an effective electric field 
within each barrier in order to satisfy the continuity equa- 
tion for the current [l5| . 

However, actual MTJ contain large amounts of disor- 
der in the electrodes, in the barriers, in the quantum well, 
and at the electrode/barrier or electrode/quantum well 
interfaces^, 0|- This disorder may represent interdif- 
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fusion at the interfaces, interface roughness, and impu- 
rities. Experiments in single MTJ have suggested that 
disorder can affect the TMR in a critical way giving rise 
to impurity-assisted tunneling 0, 0| • While the effect of 
impurities within a single barrier has been studied the- 
oretically recently \rA , its role in DMTJ remains an 
unexplored area thus far. 

In this paper we present exact model calculations of the 
spin-dependent resonant tunneling in double MTJ struc- 
tures to study for the first time the effect of magnetic and 
nonmagnetic impurities in the magnetic middle layer on 
the TMR, spin polarization, and local current as a func- 
tion of external bias. This approach conserves the conti- 
nuity of the current in contrast to Ref. ^9] . We show that 
the impurities may induce a shift of the original QWS de- 
pending on the sign of the impurity scattering potential, 
the impurity position, and the symmetry of the original 
QWS. These effects can tune selectively the spin chan- 
nels giving rise to a wide variety of novel and interesting 
spin-dependent transport phenomena, such as a dramatic 
enhancement or suppression of the TMR and the spin 
polarization, and a sign reversal of the spin polarization. 
Interestingly, the calculations reveal that even though the 
effect of the impurity on the average spin current is small 
for antisymmetric QWS, the local spin current exhibits 
strong variation. 

We employ the free-electron band model to describe 
the electron tunneling in a double MTJ with impurities 
in the metallic middle layer of width 6, whose scatter- 
ing is modeled by a spin-dependent (5-function scattering 
potential. Figure 1 shows a schematic diagram of the en- 
ergy bands for the parallel configuration under a small 
external bias for the majority and minority spin carriers, 
respectively. The geometric parameters of the DMTJ are 
chosen in such a way that the spin-dependent energies of 
the QWS, E'^, fall within the majority (minority) spin 
band and they are of antisymmetric (symmetric) charac- 
ter, respectively. The impurity potential V relative to 
the bottom of the band in the middle layer is also shown 
with green (red) line for positive (negative) sign, with the 
impurity placed at the center of the well. 
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FIG. 1: Schematic energy diagram for the parallel configu- 
ration for the majority and minority spin carriers. The en- 
ergies of the QWS denoted by fall within the majority 
(minority) band and they are antisymmetric (symmetric), re- 
spectively. The impurity 5-type potential of positive (neg- 
ative) sign relative to the bottom of the band in the well is 
denoted with red (green) line, respectively. There is no shift of 
the antisymmetric QWS, whereas the symmetric QWS shift 
towards higher (lower) values if > 0(< 0). 



Employing the WKB approximation in the barrier re- 
gion, the one-electron Green's functions for the clean 
DMTJ are solutions of the Schrodinger equation, 
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where — Vf in the metallic layers, and Vf — Ui in the 
barriers. Here, is the spin-dependent potential (or the 
bottom of the band) of the i-th metal, Ui is the potential 
barrier height in the ith barrier, and Vi{z) is the voltage 
drop within the ith layer, assumed to be nonzero only in 
the barriers. The z is the coordinate perpendicular to 
the interface, k is the in-plane wave vector of energy i?, 
and rrii is the electron effective mass in the z-th layer. 
The coefficients of the wave functions in each layer are 
determined by the boundary conditions at the interfaces. 

The one-electron Green's function in the presence of 
the impurity is determined from the Dyson equation: 

G-(p,z,p',z') ^G^{p,z,p',z') + G"{p,z,po,zo)T- x 
xG%Po,z^,p',z'), (2) 

where = ^"[1 - V Gl:^{po, zq, po, zq)]-^ is the T- 
matrix, V" — V"^-^^ 6{z — zq) — , and po,zq is the 



impurity position. Here, V" {Yimp) refers to the impu- 
rity potential relative to the bottom of the band (Fermi 
energy) of the well. Depending on the width of the quan- 
tum well and the sign of , the poles of the T-matrix 
give rise to impurity-induced shift of resonances, which 
can in turn enhance or suppress selectively the current 
density of the majority or minority spin channel for a 
given magnetization orientation. 

The local current density for spin a is given by [T9l | 

f{p -po,z) = y [f{E) - f{E + eVe^t)] X 

xD''{E,p- pa,z)dE, (3) 

where f{E) is the Fcrmi-Dirac distribution hmction and 
the transmission probability 13°' is 

D^{E, p-po,z)= (^) J2 ^U, Al^,. (4) 
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Here, V z~ (l/2)(Vz — Vz) is the antisymmetric gradi- 
ent operator, Al^,^, = (1/2)[G^"*;''-G°t';''], where G^"*;'' 
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tions, respectively. The total transmission probability 
D = + + D^'^\ where D^^) is the transmission 
probabihty in the absence of impurity, and and Z?'^' 
are the transmission probabihties proportional to the T- 
matrix and to the T-matrix squared, respectively. The 
average current j'^ can be calculated from Eq. (3) but 
with D" [E^ p — poi z) replaced with its average value 



Po,z)dp, (5) 



cN, c is the uni- 



<D^E,z) >^^J D'^{E,p 

where the number of impurities Nin 
form impurity concentration, and A'^ is the total number 
of atoms in the plane. It is important to note that the 
average current is independent of z satisfying the cur- 
rent continuity equation. In the following, we model an 
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Fe/Oxide/Fe/Oxide/Fe DMTJ, with V-^ = = 
-4.58ey and = Kj^ = V^^ = -OMeV ||. The barrier 
widths are 8 Aand 2 A, respectively, the impurity con- 
centration is c = 5%, and the impurity is placed at the 
center of the well. 

In Fig. 2 we show the TMR — {jp — jap)/ Jap 6] as 
a function of the well thickness for the perfect DMTJ 
and in the presence of non-magnetic impurities with 
Vimp = 0, -1-2, —2 eV . Here, jp (Jap) is the current den- 
sity in the parallel (antiparallel) configuration in which 
the leads are ferromagnetically (antiferromagnetically) 
aligned to the middle magnetic layer. The peak in the 
TMR at b=4.3 A for the perfect DMTJ can be dra- 
matically enhanced (suppressed) by impurities with pos- 
itive (negative) V"' which can shift the original QWS. For 
b=4.3 A, the original QWS for the majority (minority) 
spin for the parallel (antiparallel) configuration are be- 
low Ep and they are antisymmetric. On the other hand. 
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FIG. 2: TMR versus the middle layer thickness 6 in the 
presence or absence of nonmagnetic impurities for different 
values of the impurity potential, Vimp- 



the QWS for the minority (majority) spin for the paral- 
lel (antiparallel) configuration are symmetric and above 
Ep. This is due to the larger (smaller) fcj. in the well for 
the majority (minority) spin band The effect of the 
impurity on the QWS can be understood in terms of a 
simple quantum mechanical model, namely a (S-function 
impurity-potential, ^ within a potential well. Sym- 
metric QWS are shifted towards higher (lower) energies 
for positive (negative) impurity potential. In contrast, 
antisymmetric QWS are not shifted. Thus, the posi- 
tive (negative) impurity-potential shifts the QWS away 
from (closer to) the Fermi energy, suppressing (enhanc- 
ing) both jp and j\p , and hence increasing (decreasing) 
the TMR. On the other hand, the j\, and jj^p are not 
affected by the impurity. 

For b= 7.24 A, in the absence of impurity the ma- 
jority and minority spin QWS for both the parallel and 
antiparallel configurations are symmetric and below Ep. 
For Vimp = +2, eV, both the V"^ and V-^ are posi- 
tive and the QWS are shifted closer to Ep. Thus, the 
current for both parallel and antiparallel configurations 
increases and the TMR decreases. For Vimp = - 2 eV 
the QWS for the minority (majority) spin for the paral- 
lel (antiparallel) configuration are lowered in energy and 
hence the corresponding currents decrease. On the other 
hand, the QWS for the majority (minority) spin for the 
parallel (antiparallel) configuration are raised in energy 
and hence the corresponding currents increase. However, 
it turns out these current components compensate each 
other and the TMR is not altered much. 

We next investigate the effect of a magnetic impu- 
rity on the spin polarization (SP) of the tunneling elec- 
trons in DMT J. While the issue of SP for single barri- 




FIG. 3: Spin polarization for the parallel configuration as a 
function of V^^^ in the presence or absence of impurities for 
different values of V^^^ for b = 4.3 A and 7.24 A, respectively. 

ers has attracted significant interest recently, both the- 
oretically [13, nil and experimentally 0, 1^, the SP in 
DMTJ in the presence or absence of disorder has not 
been addressed theoretically to the best of our knowl- 
edge. The SP for the parallel configuration is SP = 
(jp - ip)l{il+ ipM% In Fig. 3 we show the SP for 
the parallel configuration as a function of FjJ„p for three 
values of V^^^^ for b — A.i A and 7.24 A, respectively. For 

b=4.3 A, the SP is independent of V^^^^, due to the fact 
that the QWS for the majority spin are antisymmetric 
and hence j\> is not affected by the impurity. On the 
other hand, the SP depends on V^^^ because the minor- 
ity spin QWS are symmetric. In the absence of impurity, 
jp < jp 8.nd hence the SP is negative. An impurity with 
V^ > Q decreases jp and hence increases the SP, leading 
to a sign reversal of the SP for large enough V^^^ = 2eV. 

In contrast, jp increases ii V^ < and hence the SP 
decreases. 

For b= 7.24 A, the SP increases with increasing V^l^^ 
because the symmetric QWS below Ep are shifted close 
to the Fermi energy and jp increases. Note that for 
V"J„p > —2eV the SP varies weakly with V^J^^, due to 
the fact that the QWS lie in an energy range within k^T 
around Ep. On the other hand, for V^l^^^ — and 2 eV, 
the QWS states are shifted closer to Ep {V^ > 0), and 
hence j^ increases and the SP decreases. In contrast, for 
FjI^p = -2eV the QWS states are shifted away from Ep 
{V^ < 0), and hence j^ decreases and the SP increases. 

In order to understand the effect of a non-magnetic 
impurity on the distribution of the local current in its 
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FIG. 4: Local current j''{p — pa) in the center of the first 
barrier as a function of p — po for 6 = 4.3 A and for various 
values of Vimy for the parallel configuration. The nonmagnetic 
impurity is placed at the center of the well. 

vicinity we present in Fig. 4 the local current, — po)i 
at the center of the first barrier as a function of p — po 
for h = 4.3 A and for various values of Vimp for the par- 
allel configuration. The non-magnetic impurity is placed 
at the center of the well. It is interesting to note that 
while the average current for the majority spin channel 
is not affected by the impurity due to the antisymmetric 
nature of the QWS (Fig. 2), the local current is dramati- 
cally changed compared to its corresponding value for the 



pure DMT J. The effect is larger as the impurity potential 
increases. On the other hand, for the minority spin the 
local current in the vicinity of the impurity decreases for 
Vimp = 0, +2eV and increases for Vimp = —2eV, consis- 
tent with the behavior of the average current behavior 
(Fig. 1). For distances larger than 10 A the electrons do 
not feel the impurity. 

In conclusion we have presented exact model calcula- 
tions for the effect of impurities in the well on the spin- 
dependent resonant tunneling for DMTJ. To the best 
of our knowledge these are the first calculations to ad- 
dress the effect of impurity in DMTJ. The calculations 
reveal that the spin-dependent impurity scattering po- 
tential can tune selectively the majority and minority 
spin channels, giving rise to a wide variety of interest- 
ing and unusual transport phenomena. We find that the 
impurities can lead to a dramatic enhancement or sup- 
pression or sign reversal of the TMR and the spin po- 
larization. The proposed underlying mechanism which 
explains consistently the overall behavior, is the shift of 
the original quantum well states depending on (i) the 
symmetric or asymmetric nature of these QWS and (ii) 
the impurity scattering potential. Our results for the lo- 
cal current indicate that even though the effect of the 
impurity on the average spin current is small for anti- 
symmetric QWS, the local spin current exhibits strong 
variation. 
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